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Abstract 


The spatial variability of parameters of raindrop size distribution and its derivatives 
is investigated through a field study where collocated PARSIVEL* and two- 
dimensional video disdrometers are operated at six sites in Wallops Island, Virginia 
from December 2013 to March 2014. The three-parameter exponential function is 
employed to determine the spatial variability across the study domain where the 
maximum separation distance was 2.3 km. The nugget parameter of exponential 
function is set to 0.99 and the correlation distance (do) and shape parameter (So) are 
retrieved minimizing root-mean-square error, after fitting it to the correlations of 
physical parameters. Fits were very good for almost all fifteen physical parameters. 
The retrieved do and so were about 4.5 km and 1.1, respectively, for rain rate (RR) 
when all twelve disdrometers were reporting rainfall with a rain rate threshold of 
0.1 mm h*! in one-minute observations. The do decreased noticeably when one or 
more disdrometers were required to report rain. The do was considerably different 
for a number of parameters (e.g. mass weighted diameter) but was about the same 
for the other parameters (e.g. RR) when rainfall threshold was reset to 12 dB for Ka- 
band and 18 dB for Ku-band reflectivity following the expected Global Precipitation 
Measurement mission’s space-borne radar minimum detectable signals. The 
reduction of the database through elimination of a site did not alter do as long as the 
fit was adequate. The correlations of 5-minute rain accumulations were lower when 


disdrometer observations were simulated for a rain gauge at different bucket sizes. 
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1. Introduction 


The National Aeronautics and Space Administration (NASA) Global Precipitation 
Measurement (GPM) mission aims to retrieve the three dimensional hydrometeor 
size distribution of precipitation through its Dual-frequency Precipitation Radar 
(DPR) on board the GPM core satellite (Hou et. al. 2014). The retrieval algorithm 
outputs parameters of gamma model size distribution for each range bin utilizing 
Ka- and Ku-band radar measurements (Seto et al. 2013). The lowest clutter free 
range bin near the surface at nadir incidence is approximately 5 km radius and 125 
m height. One of the key uncertainties of retrieved size distributions is its spatial 
variability within a given DPR footprint. The spatial variability is a result of 
precipitation gradient within the range volume, which may not be completely 
covered by the precipitation. This variability contributes to Non-Uniform Beam 
Filling (NUBF) which results in a higher degree of uncertainty in microwave sensor 
based precipitation estimates where the instantaneous field of view is typically 


bigger than the DPR footprint (Tokay et al. 2014a). 


Scanning radars in range height indicator mode provides the most relevant data 
source to study the spatial variability of hydrometeor size distribution in the 
vertical, while gridded radar data at 1x1 km or 2x2 km resolution is a common 
resource for studying the horizontal spatial variability. The parameters of the size 
distribution that are derived from radar measurements are based on empirical 


relationships (Bringi et al. 2004). There is also an uncertainty within the radar pixel 
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as its size increases away from radar. Disdrometers are an alternative resource for 
estimating the horizontal spatial variability but are associated with point sampling 
and it is quite costly to populate a dense network with disdrometers to sample the 
area of a satellite footprint. Therefore, there have been relatively few field studies 
conducted to determine the spatial variability of hydrometeor size distributions and 


all of these dealt with rainfall. 


As part of the NASA’s Tropical Rainfall Measurement Mission (TRMM) Ground 
Validation program, a number of impact type Joss-Waldvogel disdrometers were 
first deployed during a series of field campaigns and also deployed later at the NASA 
Wallops Flight Facility on Wallops Island, Virginia. The latter operation allowed for 
determination of the measurement accuracy of the disdrometers when six units 
were side-by-side (Tokay et al. 2005) and the spatial variability of raindrop size 
distribution (DSD) when three units were deployed across the Wallops Island strip, 
where the minimum and maximum separation distances were 0.65 and 1.7 km, 
respectively (Tokay and Bashor 2010). Later, four disdrometers were distributed 
across the same strip where the minimum and maximum separation distances were 
0.4 and 5.0 km, respectively (Schréer 2011). The three-parameter exponential 
function was tested to determine the spatial variability of DSD and integral rain 
parameters. Since four disdrometers provide six pairs of correlations across 5 km, 


the fitted exponential function was subject to noticeable error. 
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Measurement accuracy is one of the key uncertainties in determining the spatial 
variability of DSDs quantitatively. Each disdrometer type has its own shortcomings 
in measuring the DSD, and comparative field studies where different types of 
disdrometers were collocated, help quantify the uncertainties of the disdrometers 
(Krajewski et al. 2006, Thurai et al. 2011, Tokay et al. 2013, 2014b). Indeed, the 
Miriovsky et al. (2004) pioneer study was unable to determine the spatial variability 
of reflectivity within 1 km? due to uncertainties of four different types of 
disdrometers. Lee et al. (2009), on the other hand, used four Particle Occurrence 
Sensor System (POSS) to study the spatial variability of DSD in stratiform rain 
events. An S-band dual polarization radar scanned over the POSS units and the 
spatial correlation of rainfall was higher in radar than in POSS especially at 15 and 
30 km. This could partly be due to the differences between the sampling volumes of 


the instruments. 


The first comprehensive field study to quantify the spatial variability of the DSD was 
conducted in Central Spain where eight dual PARSIVEL (PARticle Size VELocity) 
disdrometers were deployed (Tapiador et al. 2010). The availability of number of 
disdrometers allowed 28 pairs of correlations where the distances range from 0.2 
km to 3.2 km. The dual units were aligned in both North-South and East-West 
directions, and one of the five events had significant differences in correlations of 
reflectivity between North-South and East-West aligned disdrometers. Jaffrain et al. 
(2011), on the other hand, analyzed 53 hours of 16 PARSIVEL observations to 


determine the spatial variability of DSD in a radar pixel (about 1 x 1 km?) in 
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Lausanne, Switzerland. They concluded that the coefficient of variations for the 
mass weighted diameter (Dmass), total concentration (Nr), and rain rate (RR) were 
high and could not be explained solely by the uncertainty of the measurement. 
Jaffrain and Berne (2012, JB12 hereafter) used the same set up of instruments to 
determine the spatial variability of Dmass, Nr, and RR using a variogram analysis. 
They reported that the variability was greater in convective rain than transitional 
and frontal rain. They also noted a decreasing variability with decreasing temporal 
resolution. Perhaps, the main issue of these studies is the measurement accuracy of 
the PARSIVEL disdrometer. The low-cost laser resulted in overestimation of large 
drops due to an inhomogeneous beam (Thurai et al. 2011, Tokay et al. 2013). The 
manufacturer upgraded the PARSIVEL disdrometer with PARSIVEL? in 2011 (Tokay 


et al. 2014b). 


In contrast to disdrometers, rain gauges are low-cost, durable, easy to maintain and 
are frequently deployed to adjust or validate radar rainfall estimates during either 
two-month long field studies or longer-term field observations. Thus, rain gauge 
networks are an excellent resource to study the spatial variability of rainfall (Habib 
and Krajewski 2002, Gebremichael and Krajewski 2004, Ciach and Krajewski 2006, 
among others). The spatial variability of rainfall is often quantified applying a three- 
parameter exponential function to the corrections between the paired gauge 
measurements. The parameters of the exponential function differs from one study 
to another due to the differences in quality of rain gauge data, inter-gauge distances, 


sample size, experiment period and location (Villarini et al. 2008). Logistics may 
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limit the inter-gauge distances, while long-term observations with dual or multiple 
gauge at each site has a clear advantage for a larger sample size and continuous 
gauge record at a given gauge site. The long-term observations also allow 
investigation of season-to-season differences in the spatial variability of rainfall. 
Tokay et al. (2014a), for instance, used five years of continuous rain gauge 
measurements in the Southern Delmarva Peninsula and showed that weather 
systems dominate spatial variability over a season. Nor’easters, for instance, 


dominate winter precipitation in the Mid-Atlantic region. 


As a comparative study, the spatial variability of rainfall was quantified through X- 
band radar and rain gauge network (Moreau et al. 2009). The correlations between 
the paired radar pixels at 1 km? resolution were significantly higher than those 
between paired gauges underneath the radar pixel up to 8 km. The differences in 
correlations are attributed to the differences in sampling volume, while the errors in 
radar rainfall estimate and time-height ambiguity also be played a role. Another 
comparative study of the spatial variability of rainfall was conducted using S-band 
dual-polarization radar and 2DVD measurements during Mid-latitude Continental 
Convective Clouds Experiment (MC3E) (Bringi et al. 2015). A good agreement was 
found between the radar and disdrometer based spatial variability for RR, median 
volume diameter, and logarithmic normalized intercept parameter with respect to 
liquid water content (Nw) in a long lasting stratiform event. The correlation 
distances were lower in a relatively short convective event than the stratiform event 


for all three parameters. 
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This study investigates the pixel-space variability of DSD and integral rain 
parameters employing a disdrometer network at NASA Wallops Flight Facility 
(WFF), Wallops Island, Virginia (37.9 degree N and 75.4 degree W). The manuscript 
is organized as follows: The disdrometer network and the database are presented in 
Section 2. Section 3 introduces the three-parameter exponential function while the 
parametric form of the DSD is given in Section 4. The probability and cumulative 
distributions of the DSD and rain parameters can be found in Section 5. Section 6 
depicts the spatial variability of DSD and rain parameters, while the sensitivity 


studies are in Section 7. The conclusions are given in the last section. 


2. Field Study 


The disdrometer network consists of six sites where each site had one two 
dimensional video disdrometer (2DVD) and one PARSIVEL? disdrometer. The sites 
were distributed across WFF where the minimum and maximum separation 
distances were 0.5 and 2.3 km, respectively (Figure 1a). One of the sites (Pad) 
hosted a variety of tipping and weighing bucket rain gauges including two pit gauges 
(Figure 1b). Another site was collocated with the Wallops Automated Surface 
Observing System (ASOS) as well as additional tipping bucket rain gauges (Figure 


1c). Table 1 summarizes the locations and distances between the sites. 
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PARSIVEL? is the third generation of the laser-optical PARSIVEL disdrometer 
(Tokay et al. 2014b). It is designed to measure the size and fall velocity of individual 
hydrometeors across its laser beam and is also present weather sensor. For rain, 
PARSIVEL outputs drop counts in a 32 x 32 size versus fall velocity matrix at 
selected time intervals ranging from 10-second to one-minute. The first two size 
bins are empty due to low signal to noise ratio and the smallest measurable 
raindrop corresponds to 0.25 mm diameter. The width of the size bins increases 
with the size of raindrops from 0.125 mm for drops to 1.2 mm, to 0.25 for drops up 
to 2.5 mm, to 0.5 mm for drops up to 5 mm, and to 1.0 mm for the larger drops. 
Therefore, there is noticeable uncertainty in determining the maximum drop size of 
the DSD. Based on 2DVD observations, the largest raindrop ever recorded is 9.7 mm 
diameter (Gatlin et al. 2015). The 25" bin corresponds to 9-10 mm and was 
considered as the largest bin for rain. For fall velocity, PARSIVEL has a range from 
0.05 to 20.8 m s“1, covering fall speeds of all types of hydrometeors. For rain, the 
26‘ bin corresponds to 8-9.6 m s!, and covers the expected range for terminal fall 
speeds of very large raindrops (Beard 1976). However, PARSIVEL? underestimates 
the fall velocity of raindrops by approximately 1 m s*! with respect to the expected 
terminal fall speed of raindrops at around 1 mm diameter (Tokay et al. 2014b). The 
underestimation in fall velocity is also evident for larger drops but the difference 
between the mean PARSIVEL? and terminal fall speed decreases with increasing size. 
The manufacturer recognizes this matter as a software error and it is expected that 
the new generation of instruments will mitigate this issue (Kurt Nemeth, OTT 


PARSIVEL, personal communication, 2015). 
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The 2DVDs are the “compact” third generation and are composing of sensor unit and 
indoor personal computer (Schonhuber et al. 2007). The third generation compact 
version has the optical components firmly mounted with no need for re-alignment 
by the user. The distance of the optical slits to the rim of the housing was also 
reduced aiming to eliminate inhomogeneous filling of the measurement area in 
windy conditions. The high-speed line-scan cameras provides better matching of 
falling hydrometeors between the two planes that are nominally 6 mm apart. The 
2DVD records the time stamp of each raindrop including its equivalent diameter, fall 
velocity, and oblateness as well as the sampling area. Common interruption of data 
occurs often due to non-meteorological items (e.g. leaves) in the sampling cross 


section, which is the main shortcoming of the 2DVD. 


The dataset for this study spanned from December 2013 through March 2014. 
While PARSIVEL? disdrometers operated nearly continuous during the experiment 
period, a number of 2DVD units failed to operate in a few rain events. The raw 
outputs of PARSIVEL* and 2DVD observations were integrated to one-minute after 
screening secondary and mismatched drops that fell outside +50% of their terminal 
fall speed (Tokay et al. 2013). Several rain/no-rain thresholds were then applied to 
the one-minute observations. All thresholds required a minimum of 10 drops. The 
minimum RR of 0.1 mm h'! resulted in 447 one-minute samples when all twelve 
disdrometers reported rainfall. Considering the minimum detectable reflectivity of 


the GPM DPR, the Ku-band reflectivity (Zxu) of 18 dB and Ka-band reflectivity (Zxa) 
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of 12 dB (Seto et al. 2013) are the other two thresholds used in this study. The 
minimum RR was 0.11 and 0.08 mm h'! for Zxa > 12 dB and 0.16 and 0.30 mm h*! for 
Zku > 18 dB for 2DVD and PARSIVEL?, respectively. The sample sizes were 445 and 
278 for Zxa > 12 dB and Zxu > 18 dB, respectively. The sample of Zxa > 12 included 


twelve one-minute spectra that were not in the RR > 0.1 mm h*! based sample. 
3. Methodology 


A three-parameter exponential function is employed to quantify the spatial 
variability of the DSD. It has been widely used to determine the spatial variability of 
rainfall through rain gauge and radar rainfall studies (Habib and Krajewski 2002, 
Gebremichael and Krajewski 2004, Ciach and Krajewski 2006, Villarini et al. 2008, 
Moreau et al. 2009, Tokay et al. 2014a). The Pearson correlation coefficient, r, is 
applied to the selected DSD and rain parameters that are derived from disdrometer 
measurements at distance d. The parametric form of exponential function is then 


expressed as 
Tai ext =" | (1) 


where ro is the correlation of a selected parameter derived from collocated 
instruments and is known as the nugget parameter and ideally should be one. 


However, collocated gauge and disdrometer observations show that the nugget 
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parameter ranges between 0.90 and 0.99. The other parameters of the exponential 
function are not sensitive to changes in ro within this range (Schroer 2011, Tokay 
and Ozturk, 2012). In this study, ro was set to 0.99. The correlation distance do, and 
the shape parameter So, are then derived through minimizing the root mean square 
error between the observed and derived correlations. As an initial guess, do values 


of 1 to 300 km and So values of 0.1 to 2.0 were given. 


The do decreases with increasing so when the correlation coefficient is high (Figure 
2). When the variability is investigated over a relatively small domain as in this 
study, the correlations of the selected DSD or rain parameters may not have a 
decreasing trend with distance. If the correlations were high, the best fit results in 
very high do and this should not be interpreted independent of other exponential 
function parameters. The parameters of exponential function are valid within the 


maximum distance of the study domain (Tokay and Oztiirk 2012). 


4. Raindrop Size Distribution 


The normalized gamma distribution function is adopted to determine the spatial 
variability of DSD parameters. The normalization was done with respect to Nr and 
liquid water content W (Tokay and Bashor 2010). The normalized intercept 
parameters with respect to total concentration, Nr’ and liquid water content, Nw are 


expressed as 
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where pw is the density of water. 


(2) 


(3) 


Dmass is related to the slope, A and shape 


parameter, m of the complete gamma distribution as; 


(4) 


The normalized intercept parameters can then be calculated from observed spectra 


as well. The corresponding normalized gamma fitted distributions are expressed as: 


ND) = Ni, (m(— 


mass 


N(D) = Nyf(m)(— 


mass 


where fi(m) and f2(m) are given as; 


(4 + m)™*! 


n= Tim+1) 


)™ exp[-(4 + m) 


)™ exp[-(4 + m) 


| (5) 


mass 


| (6) 


mass 


(7) 
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(8) 


To extract the shape parameter, m(N, ), and m(N,), the rain rates that are calculated 
from observed and fitted gamma distributions (e.g. equations 5 and 6) are 
minimized. The formulations presented above follows the complete gamma 
function where minimum and maximum drop size is assumed zero and infinity, 
respectively. In reality, there is a minimum and maximum drop size (Dmax) in a 
population of drops and the incomplete gamma function is more appropriate 
especially if the size spectra do not contain the large drops (e.g. narrow 


distribution), but this is the beyond the scope of this study. 


The GPM DPR algorithm has adopted the normalized gamma function as in equation 
6 (Seto et al. 2013). Since the DPR algorithm employs dual frequency reflectivity 
measurements to determine the three parameters of gamma function, there is an 
interest in finding relations between the derived parameters. Williams et al. (2014) 
suggested a power law relation between the standard deviation of Dmass (Omass) and 
Dmass. Both variables are directly calculated from observed DSD and their ratio is a 


sole function of shape parameter, m(Omass) following complete gamma function. 


Omass _ 1 (9) 
Dynass (4+ m)?> 
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5. Probability and Cumulative Distributions 


This study investigates the spatial variability of seven DSD and eight integral rain 
parameters. The DSD parameters include Dmass, Dmax, Nr’, Nw, m(Nr’), m(Nw), and 
m(Omass), While the rain parameters were W, RR, dual polarization parameters of 
horizontal reflectivity (Zu), and differential reflectivity (Zar), reflectivity at W-band 
(Zw), Zku, Zka, and dual frequency ratio (DFR = Zxu/Zxa). The Ka-, Ku-, and W-band 
reflectivities are calculated using Mie scattering for spherical particles, while dual- 
polarization parameters are calculated for S-band radar following Tokay et al. 
(2002). Table 2 presents the mean, standard deviation, median, and 5t and 95 
percentage of 15 physical parameters. Both Nr’ and Nw have ranges of several order 
of magnitude and, hence logarithmic values of these two parameters are used in 


constructing probability and cumulative distributions. 


Among many factors listed in the introduction, knowledge of the characteristics of 
the DSD and rainfall is essential for quantifying the spatial variability. Observed 
DSD lacking large drops but with abundant small drops will incur different spatial 
variability than a DSD with numerous large drops but fewer small drops. The 
spatial variability of rainfall, on the other hand, differs in the presence and absence 
of heavy rain. The probability and cumulative distributions of the DSD and rain 


parameters provide an insight on the characteristics of DSD and rainfall, and these 
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distributions should be included if a similar study is conducted with a different 


dataset. 


Figure 3 presents the probability and cumulative distributions of fifteen DSD and 
integral rain parameters based on 2DVD and PARSIVEL? measurements, when all 
twelve units were reporting rainfall. A very good agreement in midsize range 
(diameter 1-3 mm) of 2DVD and PARSIVEL? size spectra reflects an excellent 
agreement in probability and cumulative distributions of Dmass between the two 
disdrometers (Figure 2a). Dmass mostly resided between 0.8 and 1.5 mm, peaking at 
1 mm (Table 2). The binning of drop counts results in quantization error in size 
measurements in PARSIVEL?, which is quite significant for Dmax since the larger bins 
have wider widths. The probability distribution of Dmaxin PARSIVEL? therefore had 
multiple modes, while the probability distribution of Dmaxin 2DVD was a unimodal 
(Figure 3b). The cumulative distributions of Dmax agreed well between the two 


types of disdrometers and the median Dmax was around 1.9 mm for both (Table 2). 


PARSIVEL? is relatively more sensitive to small drops less than 0.5 mm diameter 
than the 2DVD (Tokay et al. 2013, 2014b). Nr’ is very sensitive to the number of 
small drops since their concentrations are much higher than larger drops. The 
probability distribution of log(Nr’) of 2DVD and PARSIVEL? show a wide range, but 
the PARSIVEL* based distribution was shifted to the larger concentrations, 
consistent with its sensitivity to small drops (Figure 3c). Nw is more sensitive to 


midsize drops and therefore an agreement is expected between the 2DVD and 
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PARSIVEL? distributions. A general agreement is evident between the distributions 
of log(Nw) with a median 3.55, but the PARSIVEL? also had low percentages of 


log(Nw) less than 2.5 (Figure 3d). 


The shape parameter of the gamma distribution is very sensitive to the method of 
derivation and has a wide range but was mostly between 0 and 14 with mean and 
median around 5 (Table 2). While the agreement between 2DVD and PARSIVEL2 
derived distributions of shape parameters were reasonable for m(Nr’) m(Nw), and 
m(Omass), the probability distribution of m(Nr’) was shifted toward larger values in 
PARSIVEL? and m(Nw) and m(Omass) had slightly larger values in 2DVD (Figures 3e- 
g). The differences in probability and cumulative distributions are more 


pronounced in m(Nr’) reflecting the differences in distributions of log(Nr’). 


An excellent agreement was evident in the distributions of W, and RR between the 
two different types of disdrometers (Figures 3h-i). Both W and RR are primarily 
sensitive to midsize drops, while the upper end of the small drops (0.8-1.0 mm) and 
the lower end of the large drops (3.0-4.0 mm) contribute to W and RR significantly 
in the presence of light and heavy rainfall, respectively. The mean and median RR is 
1.1 and 0.9 mm h'!, respectively (Table 2), indicating the dominance of light rain. 
Good agreement was also evident in distributions of Zu, Zar, Zw, Zku, Zka, and DFR 
between the two types of disdrometers (Figures 3j-o). The median Zu, Zxu and Zxa 
were in the 22-23 dB range, while 95" percentiles were in the 31-32 dB range 


(Table 2) in the presence of widespread frontal rainfall. 
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6. Spatial Variability 


The correlations of Dmass and Dmax are higher for 2DVD than for PARSIVEL? at a given 
distance (Figures 4a-b). PARSIVEL? based correlations were 0.55-0.7 for five paired 
observations resulting in a higher RMSE. The quantization error due to binning 
contributed to the low correlations in PARSIVEL? based Dmax. The correlations of 
Nr’ and Nw were higher than 0.9 for all distances for the 2DVD and decreased 
gradually with distance (Figures 4c-d). PARSIVEL? based correlations have a wider 
range than 0.1 at a given distance but remained higher than 0.84 and the RMSE was 
quite low for fitted exponential function. The correlations of shape parameters 
decreased with distance from 0.8 to 0.5 for m(Nr’) and from 0.8 to 0.4 for m(Nw) 
and m(Omass) for both 2DVD and PARSIVEL? (Figures 4e-g). The range of 
correlations was less than 0.1 at a given distance for both disdrometers resulting in 
low RMSE. The sample size was 74% and 77% of the database for m(Nr’) and 
m(Nw), respectively due to the fact that the shape parameter was outside the 
expected range of -4 to 20 for a number of size spectrum when normalized gamma 
function was fitted minimizing rain rate. Similarly, the sample size was 76% of the 
database for m(Omass) due to presence of shape parameters larger than 20. The 
large shape parameters results from narrow size spectrum in the absence of large 
drops. The shape of size spectrum exhibits more than a single peak with a plateau 


in midsize regime for collisional break-up dominated DSD (D’Adderio et al. 2015). 
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The gamma distribution is not the best mathematical fit in these conditions, but this 


is beyond scope of this paper. 


The correlations of RR and W decreased with distance from 0.9 to 0.55 and from 
0.92 to 0.64, respectively (Figures 4h-i). The spread in correlation between the 
same and different types of disdrometers was much less than 0.1 at a given distance 
resulting in very low RMSE. Thus, the exponential fits were nearly identical for both 
types of disdrometers. The correlations of reflectivities at different wavelengths 
decreased with distance but there were noticeable differences and similarities 
among them. The correlations of Zy and Zxu were about the same at a given distance 
for a given disdrometer while the spread in PARSIVEL? derived correlations of Zxa 
were less resulting in lower RMSE (Figures 4j, 4m-n). The spread in correlations of 
Zw were much less (less than 0.05) at a given distance for both disdrometer types 
resulting in an excellent fit (Figure 41). For Zar, there were significantly lower 
correlations (0.35-0.55) for five PARSIVEL? pairs resulting in high RMSE (0.14) 
(Figure 4k). For 2DVD, Zar decreased with distance and had a spread of 0.1 at a given 
distance resulting in a reasonable fit with low RMSE. DFR, which is another 
reflectivity ratio, had lower correlations even at short distances and the spread in 
correlations were as high as 0.3 at a given distance in PARSIVEL? resulting in 
relatively poor fit with high RMSE (0.12) (Figure 40). For 2DVD, the fit was better 


but the spread in correlation was 0.2 at a given distance. 
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The correlation distances were about the same for RR and W but higher in 2DVD 
than in PARSIVEL? for the rest of rain and DSD parameters (Figure 5a). For Nr’ and 
Nw, do was at around upper limit of initial guess due to high correlations for all 
distances. The do was 4.2 and 4.5 km in RR in PARSIVEL? and 2DVD, respectively. 
These correlation distances correspond to high spatial variability following gauge- 
based studies (Ciach and Krajewski 2006). However, the correlations remained 
high (> 0.7) in this study (Figure 4i). The gauge studies were conducted for 5- 
minute or longer integration periods while high temporal scale in this study results 
in higher variability even in the presence of lighter rain. The shape parameter, on 
the other hand, remained between 0.4 and 1.0 for most of the fits (Figure 5b). The 
RMSE was less than 0.08 for almost all parameters showing the goodness of fit 
(Figure 5c). It should be emphasized that three-parameter exponential function is 
applied for the first time for a number of DSD and rain parameters following 
feasibility study (Schréer 2011). Therefore, there is no direct comparison available 


for these parameters of exponential function with any other studies. 


7. Sensitivity Studies 


a) Disdrometer network availability 


It is not uncommon that one or more disdrometers fail to operate throughout a field 


campaign. For spatial variability, it is important to assess the goodness of fit of the 


three-parameter exponential function when one of the sites is not available. In that 
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regard, five out of six sites were employed to quantify the sensitivity of parameters 
of exponential function to Dmass and RR. The elimination of a site reduced paired 
correlations from 15 to 10 and if site 5 or site 6 in Table 1 were not functioning, the 


maximum separation distance would be reduced to 1.84 km. 


Fits were very good for 2DVD derived Dmass and RR regardless of the choice of 
elimination of paired correlations. The do had a very narrow range of 4.2-4.7 km for 
RR but remained mostly within 11-21 km for Dmass (Figure 6a). For PARSIVEL2, fits 
were good for RR but quite poor for Dmass repeating previous finding in Figure 4a. 
The do had a wide range in PARSIVEL? derived Dmass but was not sensitive to the 
elimination of any paired correlation in RR. The so was less than 1.0 and highly 
variable due to elimination of paired correlation in Dmass but remained mostly above 
1.0 with almost no sensitivity to the elimination of paired correlation in RR (Figure 
6b). The RMSE was 0.03 or less for 2DVD for both Dmass and RR but was higher than 


0.1 for four trials for PARSIVEL? derived Dmass (Figure 6c). 


b) Rain/No-rain threshold 


The rain/no-rain threshold results in differences in the sample size and the 
distribution of the DSD and rain parameters. The higher sensitivity DPR footprint 
where Zxa is larger than 12 dB is expected to have less spatial variability with higher 
do than the lower sensitivity DPR footprint where Zxu is larger than 18 dB. As shown 


in the gauge based study of Tokay and Ozturk (2012), the higher thresholds 
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eliminate the light rain or low reflectivity samples and results in lower correlations, 


especially if there is no precondition where both disdrometers must report rainfall. 


Based on 2DVD observations, do of Dmass was the lowest and the highest for Zxu and 
RR thresholds, respectively (Figure 7a). The do of RR and W was the same for RR 
and Zxa thresholds and slightly lower for Zxu threshold. The do of Za, Zku, and Zu 
were about the same for Zxa and Zxy thresholds, and were lower than RR threshold. 
The do of DFR was the same for all three thresholds, while the other DSD and rain 
parameters showed different trends. The shape parameter of exponential function 
showed little or no variability when do was about the same for given thresholds 
(Figure 7b). The RMSE did not show significant differences between the different 
thresholds for most of the DSD and rain parameters except for Zar and DFR where 


RMSE was noticeably higher for Zxu threshold (Figure 7c). 


c) Rain coverage 


Partial coverage of satellite footprint or radar pixel by rain is one of the sources of 
NUBF and contributes significantly to the spatial variability. To quantify the spatial 
variability in the presence of partial coverage, PARSIVEL? observations were 
reprocessed when all units were reporting rainfall and when at least one unit 
reporting rainfall. The latter is more commonly observed in nature. The sample 
size is 4,645 one-minute size distributions for the latter condition, which makes it 


more than 2.6 times of the sample size of the former condition. The difference is 
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mainly due to the rain intermittence. The sample size when all six PARSIVEL2 
reported rainfall was 4 times larger than that when all twelve disdrometers 


reported rainfall. 


The correlation distance was distinctly lower for almost all physical parameters 
when one or more PARSIVEL? reported rainfall (Figure 8a). It should also be noted 
that do was different from earlier findings when all six PARSIVEL? were reporting 
rainfall. The bigger sample included more uniform rainfall across the study domain 
and do of RR was 14 km. The so was higher when one or more PARSIVEL? reported 
rainfall (Figure 8b). Fits were good for both conditions for physical parameters 
except for normalized intercept and shape parameter of the gamma distribution 


where fit was relatively poor with RMSE of 0.1 or higher (Figure 8c). 


d) Rain gauge simulation 


Tipping bucket rain gauges are widely used in precipitation studies but suffer from 
significant sampling errors over short integration periods depending on the bucket 
size (Habib et al 2001). A bucket size of 0.01 inch (0.254 mm) is used in ASOS 
network but the gauge manufacturers also provide gauges with bucket resolutions 
of 0.1 mm and 0.2 mm. The gauge based spatial variability of rainfall studies 


employs 5-minute or longer integration periods to mitigate the sampling errors. 
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PARSIVEL? and 2DVD one-minute rain rate time series were employed to simulate 
the time of tip for 0.1 mm, 0.2 mm, and 0.254 mm bucket resolutions. The 
simulation study ignored any disdrometer malfunctions throughout the experiment 
period. The 5-minute rainfall was then calculated from one-minute disdrometer 
and three different simulated gauge datasets. The nugget parameter of 0.99 was 
considered in all simulations. The correlations decrease noticeably with the size of 
the bucket at a given distance in both PARSIVEL? and 2DVD based simulations 
(Figure 9a-b). The decrease in correlations was gradual for disdrometer but quite 
sharp for coarser relation bucket size based simulations from nugget to short 
separation distances. At 2 km separation distance, the differences in correlations 
were approximately 0.2 between disdrometer and 0.254 mm bucket simulations of 


2DVD and PARSIVEL?. 


8. Conclusions 


The pixel-scale variability of seven DSD and eight integral rain parameters was 
investigated through a unique set of disdrometer observations where collocated 
2DVD and PARSIVEL? were operated at six sites across the main base of NASA/WFF. 
A three-parameter exponential function was employed to quantify the spatial 
variability. The ro, or nugget parameter, was fixed to 0.99 and do and so were 
retrieved minimizing RMSE following fitting an exponential function to the 
correlations that were derived from 2DVD or PARSIVEL? observations. The 


correlations were calculated during periods when all twelve disdrometers reported 
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rainfall and a minimum RR threshold of 0.1 mm h’! was satisfied with a minimum of 


10 drops occurring in one-minute. 


The there parameter exponential function is a simple mathematical form that 
expresses the observed correlations well. It converged rapidly for almost all 
physical parameters where RMSE was mostly less than 0.08. Very good agreement 
between the 2DVD and PARSIVEL? derived correlations boosted our confidence in 
the quality of the observations. The 2DVD outperformed the PARSIVEL? for 
sampling of Dmax, while the reverse was true for Nr. This is likely due to the 
difference in the sensitivity of the respective disdrometers to the small and large 


drop end of the size spectrum. 


Given the confidence in measurement quality and fitting method, do and so should 
mainly be sensitive to the characteristics of the DSD. The exclusion of a site did not 
change the distribution of rainfall resulting in insignificant changes in do and so of 
RR. The retrieved parameters of exponential function fits quite different in 
PARSIVEL? derived Dmass when one site was removed but this is attributed to the 
poor fitting. The do and so were sensitive to the different rain/no-rain thresholds. 
They were also sensitive to the differences in conditions where one or more and all 
six PARSIVEL? were reporting rainfall. The differences in the sample of 
observations play a key role for differences in do and so. The disdrometers have a 
clear advantage to the rain gauges when spatial variability was studied for short 


time intervals. The larger the bucket size the lower the correlations when the 
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disdrometer observations were used to simulate rain gauge observations at 


different bucket size. 


This is the first study where the three-parameter exponential function has been 
used for estimating the variability of 15 physical parameters. Previously, JB12 used 
the exponential function for RR. Their study was more robust in terms of the 
number of available instruments and a longer period of observation. The 
correlations decreased more rapidly for the first 500 m reaching 0.8 in JB12 study. 
This study did not have any paired correlation in this short distance regime and the 
correlation of RR was about 0.9 at 500 m. The differences in correlations are 
attributed to the nature of rain. JB12 study included higher rainfall intensities than 


this study. 


The logistics often dominate the experiment set up and therefore the minimum and 
maximum distances in the study domain. The sample size and type of rainfall play a 
crucial role in quantifying the spatial variability. WFF is ideal site since rainfall from 
multiple weather systems falls in all around year. In this study, light rain was 
persistent throughout field study where mean RR was 1.1 mm h*! and the 95th 
percentile of Dmax was 3.1 mm (2DVD). For future studies, the larger domain where 
the maximum disdrometer distance exceeds the DPR footprint and the combination 
of light and heavy rain is desirable for quantification of spatial variability of DSD. 


The limited but unique data set from Mid-latitude Continental Convective Clouds 
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experiment from Northern Oklahoma may prove to be a great resource for such a 


follow-up study. 
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Figure Captions: 


Figure 1. Google Map of six field sites at the NASA Wallops Flight Facility (top). 
Picture of 2DVD, PARSIVEL* among other precipitation measuring instruments at 
the Pad (middle) and at the Automated Surface Observing System site (bottom). 
Please note that not all the instruments collecting data at the Pad during the 


experiment period. 


Figure 2. Dependence of the correlation distance to correlation coefficient at a given 
distance and shape parameter when nugget parameter is set to 0.99. Four different 


shape parameters and two different distances are used. 


Figure 3. Probability and cumulative distributions of (a) mass weighted diameter, 
(b) maximum diameter, (c) logarithmic normalized intercept parameter, Nr’, (d) 
logarithmic normalized intercept parameter, Nw, (e) shape parameter with respect 
to Nr’, (f) shape parameter with respect to Nw,, (g) shape parameter with respect to 
Omass, (h) liquid water content, (i) rain rate, (j) horizontal reflectivity, (k) differential 
reflectivity, (1) reflectivity at W-band, (m) reflectivity at Ku-band, (n) reflectivity at 
Ka-band, (0) dual frequency ratio. Distributions of these physical parameters are 


derived from 2DVD and PARSIVEL? observations. 


Figure 4. Spatial variability of (a) mass weighted diameter, (b) maximum diameter, 


(c) logarithmic normalized intercept parameter, Nr’, (d) logarithmic normalized 
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intercept parameter, Nw, (e) shape parameter with respect to Nr’, (f) shape 
parameter with respect to Nw,, (g) shape parameter with respect to Omass, (h) liquid 
water content, (i) rain rate, (j) horizontal reflectivity, (k) differential reflectivity, (1) 
reflectivity at W-band, (m) reflectivity at Ku-band, (n) reflectivity at Ka-band, (0) 
dual frequency ratio, all derived from 2DVD (blue dots) and PARSIVEL? (red stars) 
observations. The parameters of the three-parameter exponential function 


including root-mean-square error and the sample size are also given. 


Figure 5. (a) Correlation distance and (b) shape parameter of the three-parameter 
exponential function and (c) root-mean-square error for fifteen physical parameters 
based on 2DVD (blue dots) and PARSIVEL? (red stars) observations. Several 
correlation distances were higher than y-axis range and are marked with their 


values. 


Figure 6. Sensitivity of (a) the correlation distance, (b) shape parameter of the 
three-parameter exponential function and of (c) root-mean-square error of mass 
weighted diameter (left) and rain rate (right) to the elimination of a site (e.g. site1) 
during the experiment. Several correlation distances were higher than y-axis range 


and are marked with their values. 


Figure 7. Sensitivity of (a) the correlation distance, (b) shape parameter of the 
three-parameter exponential function and of (c) root-mean-square error of fifteen 


physical parameters to the rainfall threshold following 2DVD observations. The 


34 


750 


751 


752 


753 


754 


755 


756 


757 


758 


759 


760 


761 


762 


763 


764 


rainfall thresholds of Zxa > 12 dB (green dot), Zku > 18 dB (black dot), and RR > 0.1 m 
h-! (red dot) are considered. Several correlation distances were higher than y-axis 


range and are marked with their values. 


Figure 8. Sensitivity of (a) the correlation distance, (b) shape parameter of the 
three-parameter exponential function and of (c) root-mean-square error of fifteen 
physical parameters to rain coverage. One or more (orange stars) as well as all six 
(red stars) PARSIVEL? reporting rainfall was considered. Several correlation 


distances were higher than y-axis range and are marked with their values. 


Figure 9. Spatial variability of 5-minute rainfall derived from (a) 2DVD and (b) 


PARSIVEL? observations (red), simulated gauge at 0.1 mm bucket (blue), simulated 


gauge at 0.2 mm bucket (green), and simulated gauge at 0.254 mm bucket (black). 
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765 Table 1. NASA Wallops Flight Facility Precipitation Data Acquisition Network. The 
766 distances between the six sites are given in km. 
767 
Site # Site 1 Site 2 Site 3 Site 4 Site 5 Site 6 
Site Pad ASOS Balloon | Building Visitor Water 
name Launch A41 Center | Treatment 
Facility Plant 
Coord | 37.944°N, | 37.944°N, | 37.938°N, | 37.934°N, | 37.929°N, | 37.937°N, 
inates | 75.464°W | 75.481°W | 75.456°W | 75.471°W | 75.473°W | 75.466°W 
Site 1 0.51 0.60 1.25 1.34 1.41 
Site 2 1.06 1.81 0.88 1.53 
Site 3 1.84 1.75 1.80 
Site 4 0.99 1.52 
Site 5 2.31 
768 
769 
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770 Table 2: Mean, standard deviation, median, and 5th and 95th percentiles of seven 
771 DSD and eight integral parameters that were derived from 2DVD and PARSIVEL? 
772 (Par?) observations. 
773 
DSD and mean standard median 5% percentile | 95* percentile 
integral rain deviation 
parameters 
2DVD | Par | 2DVD Par2 2DVD | Par? | 2DVD | Par? | 2DVD Par2 
Dmmass (mm) 1.11 1.09 0.30 0.32 1.05 1.01 0.72 0.72 1.64 1.63 
Dmax (mm) 2.03 1.96 0.58 0.63 1.92 1.93 1.28 1:72 3.07 3.35 
Nr (mmm) | 302 951 394 2,271 184 319 47 64 1,065 4,274 
Nw(m3mm*) | 6,524 | 7,313 | 9,294 | 10,915 | 3,468 | 3,588] 711 53 26,787 | 32,706 
m(Nr’) 4.85 5.95 2.92 3.87 4.40 5.10 1.10 1.40 10.30 13.65 
m(Nw) 5.16 5.00 3.38 4.01 4.50 4.00 0.80 0.10 11.40 12.90 
m(Omass) 5.75 5.41 3.54 4.18 5.23 4.60 0.90 | -0.01 | 12.14 13.45 
RR (mmh‘) 1.10 1.15 0.94 1.00 0.87 0.91 0.25 0.27 2.94 2.92 
W (gm) 0.071 | 0.077 | 0.054 | 0.056 | 0.059 | 0.065 | 0.019 | 0.021 | 0.188 0.191 
Zn (dB) 26.4 27.7 18.0 33.4 22.3 22.1 15.0 15.0 31.2 31.8 
Zar (dB) 0.376 | 0.394 | 0.291 | 0.427 | 0.287 | 0.273 | 0.119 | 0.112 | 0.839 0.938 
Zw(dB) 16.9 17.3 3.4 3.2 15.8 16.4 11.3 11.8 20.8 21.1 
Zxku (dB) 26.8 Zick 18.5 22.1 21.9 21.7 14.8 14.7 chile Bas) 
Zka (dB) 26.3 26.2 5.8 6.3 23.4 23.2 15.2 15.5 31.4 Bly 
DFR -0.99 | -0.86 | 1.19 1.11 -1.27 | -1.04 | -2.04 | -1.86 1.26 0.78 
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Figure 1. Google Map of six field sites at the NASA Wallops Flight Facility (top). 
Picture of 2DVD, PARSIVEL* among other precipitation measuring instruments at 
the Pad (middle) and at the Automated Surface Observing System site (bottom). 
Please note that not all the instruments collecting data at the Pad during the 
experiment period. 
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789 distance and shape parameter when nugget parameter is set to 0.99. Four different 
790 shape parameters and two different distances are used. 
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Figure 3. Probability and cumulative distributions of (a) mass weighted diameter, 
(b) maximum diameter, (c) logarithmic normalized intercept parameter, Nr’, (d) 
logarithmic normalized intercept parameter, Nw, (e) shape parameter with respect 
to Nr’, (f) shape parameter with respect to Nw, (g) shape parameter with respect to 
Omass, (h) liquid water content, (i) rain rate, (j) horizontal reflectivity, (k) differential 
reflectivity, (1) reflectivity at W-band, (m) reflectivity at Ku-band, (n) reflectivity at 
Ka-band, (0) dual frequency ratio. Distributions of these physical parameters are 
derived from 2DVD and PARSIVEL? observations. 
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Figure 4. Spatial variability of (a) mass weighted diameter, (b) maximum diameter, 
(c) logarithmic normalized intercept parameter, Nr’, (d) logarithmic normalized 
intercept parameter, Nw, (e) shape parameter with respect to Nr’, (f) shape 
parameter with respect to Nw,, (g) shape parameter with respect to Smass, (h) liquid 
water content, (i) rain rate, (j) horizontal reflectivity, (k) differential reflectivity, (1) 
reflectivity at W-band, (m) reflectivity at Ku-band, (n) reflectivity at Ka-band, (0) 
dual frequency ratio, all derived from 2DVD (blue dots) and PARSIVEL? (red stars) 
observations. The parameters of the three-parameter exponential function 
including root-mean-square error and the sample size are also given. 
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Figure 5. (a) Correlation distance and (b) shape parameter of the three-parameter 
exponential function and (c) root-mean-square error for fifteen physical parameters 
based on 2DVD (blue dots) and PARSIVEL? (red stars) observations. Several 
correlation distances were higher than y-axis range and are marked with their 
values. 
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Figure 6. Sensitivity of (a) the correlation distance, (b) shape parameter of the 
three-parameter exponential function and of (c) root-mean-square error of mass 
weighted diameter (left) and rain rate (right) to the elimination of a site (e.g. site1) 
during the experiment. Several correlation distances were higher than y-axis range 
and are marked with their values. 
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837 Figure 7. Sensitivity of (a) the correlation distance, (b) shape parameter of the 
838 three-parameter exponential function and of (c) root-mean-square error of fifteen 
839 physical parameters to the rainfall threshold following 2DVD observations. The 
840 rainfall thresholds of Zxa > 12 dB (green dot), Zxu > 18 dB (black dot), and RR> 0.1m 
841 hr! (red dot) are considered. Several correlation distances were higher than y-axis 
842 range and are marked with their values. 
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848 Figure 8. Sensitivity of (a) the correlation distance, (b) shape parameter of the 
849 three-parameter exponential function and of (c) root-mean-square error of fifteen 
850 physical parameters to rain coverage. One or more (orange stars) as well as all six 
851 (red stars) PARSIVEL? reporting rainfall was considered. Several correlation 
852 distances were higher than y-axis range and are marked with their values. 
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859 ‘Figure 9. Spatial variability of 5-minute rainfall derived from (a) 2DVD and (b) 
860  PARSIVEL? observations (red), simulated gauge at 0.1 mm bucket (blue), simulated 
861 gauge at 0.2 mm bucket (green), and simulated gauge at 0.254 mm bucket (black). 
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